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Abstract 

We propose a scenario in which a simple power-like primary spectrum for protons with 
sources at cosmological distances leads to a quantitative description of all the details of 
the observed cosmic ray spectrum for energies from 10 17 eV to 10 21 eV. As usual, the 
ultrahigh energy protons with energies above Eqzk ~ 4 • 10 19 eV loose a large fraction of 
their energies by the photoproduction of pions on the cosmic microwave background, which 
finally decay mainly into neutrinos. In our scenario, these so-called cosmogenic neutrinos 
interact with nucleons in the atmosphere through Standard Model electroweak instanton- 
induced processes and produce air showers which are hardly distinguishable from ordinary 
hadron- initiated air showers. In this way, they give rise to a second contribution to the 
observed cosmic ray spectrum - in addition to the one from above mentioned protons - 
which reaches beyond Eqzk- Since the whole observed spectrum is uniquely determined by 
a single primary injection spectrum, no fine tuning is needed to fix the ratio of the spectra 
below and above Eqzk- The statistical analysis shows an excellent goodness of this scenario. 
Possible tests of it range from observations at cosmic ray facilities and neutrino telescopes 
to searches for QCD instanton-induced processes at HERA. 



*On leave from Institute for Theoretical Physics, Eotvos University, Budapest, Hungary. 



1 Introduction 



The puzzle of ultrahigh energy cosmic rays (UHECRs) is about 40 years old. About twenty 
mysterious events were observed above 10 20 eV by five different air shower observatories (AGA- 
SA 0, Fly's Eye 0, Haverah Park |Sj, HiRes [t], and Yakutsk 0; for reviews, see Ref. ©). 
Though some small-angle clustering in the arrival direction of the UHECRs is observed, the 
overall event distribution is isotropic. This indicates that they originate from several sources. No 
source is known, however, within a distance of 50 Mpc. This is rather peculiar, since 50 Mpc is 
the characteristic distance ultrahigh energy nucleons travel before they loose a large fraction of 
their energy. A sharp drop around the Greisen-Zatsepin-Kuzmin (GZK) cutoff Sgzk~ 4- 10 19 eV 
is therefore predicted in the cosmic ray spectrum [Jj. The available data show no such drop 1 . 

The reason for the expected drop is a well established elementary process. Above Eqzk, protons 
produce pions through the interaction with photons from the 2.7 K cosmic microwave background 
(CMB). The produced pions decay, resulting in the so-called cosmogenic neutrinos 0. The 
attenuation length of protons above the GZK cutoff is about 50 Mpc. The basic question is: if 
the sources of ultrahigh energy cosmic rays are indeed at cosmological distances, how could they 
reach us with energies above 10 20 eV? No conventional explanation is known for this question. 

At the relevant energies among the known particles only neutrinos can propagate without signif- 
icant energy loss from cosmological distances to us. It is this fact which led, on the one hand, to 
scenarios invoking hypothetical - beyond the Standard Model - strong interactions of ultrahigh 
energy cosmic neutrinos jH] and, on the other hand, to the Z-burst scenario fUl El E] • 

In the latter, ultrahigh energy cosmic neutrinos (UHECz/s) produce Z-bosons through annihilation 
with the relic neutrino background from the big bang. On earth, we observe the air showers 
initiated by the protons and photons from the hadronic decays of these Z-bosons. Though the 
required ultrahigh energy cosmic neutrino flux is smaller than present upper bounds, it is not 
easy to conceive a production mechanism yielding a sufficiently large one. In the near future, the 
neutrino telescopes AMANDA pQ| and RICE [T3], as well as the Pierre Auger Observatory 
for extensive air showers, can directly confirm or exclude this scenario. 

Scenarios based on strongly interacting neutrinos use the fact that the observed cosmic ray flux 
above Eqzk can be fairly well described by the predicted [THJ HI] cosmogenic neutrino flux. In 
these scenarios, neutrinos with energies above ~ 10 20 eV originating from the GZK process are 
assumed to have a large cross-section for the scattering off nucleons and to initiate extensive air 
showers high up in the atmosphere, like hadrons. This is usually ensured by new types of TeV-scale 
interactions beyond the Standard Model, such as arising through gluonic bound state leptons [18J, 
TeV-scale grand unification with leptoquarks or Kaluza-Klein modes from extra compactified 
dimensions [2U| (see, however, Ref. |2T]): for earlier and further proposals, see Refs. |22J and |2~3j . 
respectively. Until now, none of these ideas have direct experimental verification. 

In this Letter, we propose another strongly interacting neutrino scenario to solve the GZK prob- 
lem, which - in contrast to previous proposals - is based entirely on the Standard Model of particle 
physics. It exploits non-perturbative electroweak instanton-induced processes for the interaction 
of cosmogenic neutrinos with nucleons in the atmosphere, which may have a sizeable cross- section 
above a threshold energy E t h = 0((Ati M w / 'a w ) 2 ) / '(2m p ) = C(10 18 ) eV, where M w denotes the 

1 There is an ongoing debate whether the excess of events above Eqzk is significant and whether the data from 
different collaborations are mutually consistent [S]. We will comment on this point below. 
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W-boson mass and aw the electroweak fine structure constant |23J I2H1 12E] • F° r the first time in 
the literature, we present a detailed statistical analysis of the agreement between observations 
and predictions from strongly interacting neutrino scenarios. 

Our scenario can be summarized as follows. We assume a standard power-like primary spectrum 
for protons injected from sources at cosmological distances, which extends beyond the GZK cutoff. 
After propagation through the CMB, the protons - arriving at earth mostly with energies below 
Eqzk - will be one component of the observed cosmic ray spectrum. The spectrum of the produced 
cosmogenic neutrinos is entirely determined by the proton injection spectrum and can therefore 
be determined precisely, including all known effects. The cosmogenic neutrinos travel unaffected 
through the CMB. However, for energies above ~ 10 19 eV, they have a large cross-section for 
interactions with nuclei in the atmosphere due to electroweak instanton-induced processes. They 
give rise, therefore, to a second, predictable component of the observed cosmic ray spectrum, 
which dominates above the GZK cutoff over the first, proton- initiated component. Our proposal 
leads to an explanation of the observed cosmic ray spectrum simultaneously above and below the 
GZK cutoff, without the need to fix the ratio of the fluxes below and above E GZ k by hand 2 , as it 
is necessary in most alternative proposals. The goodness of the scenario is studied by statistical 
methods and an excellent agreement is seen between the predictions and the observations. 

Our analysis proceeds in three steps, which are performed in Sects. 121 El and |U i) First, we 
study the consequences of a power-like proton injection spectrum. We determine the resulting 
proton and neutrino fluxes on earth, taking into account the appropriate types of energy losses. 
ii) In the second step, we calculate the spectrum of cosmogenic neutrino-initiated electroweak 
instanton-induced air showers, in) The third step consists in the comparison of the observed 
UHECR spectrum with the prediction arising from an inclusion of instanton-induced processes. 
Based on the goodness of the scenario, we determine the confidence region in the parameter space 
of our scenario. Finally, we summarize our result and present our conclusions in Sect. |3J 

2 Proton and cosmogenic neutrino fluxes 

We start with a power-like injection spectrum per co-moving volume of protons with energy Ei, 
spectral index a, and redshift (z) evolution index m, 

j p = j E~ a (1 + z) m 6»(£ max - Ei) 9(z - z mhl ) 9(z max - z) . (1) 

Here, jo is a normalization factor, E max is the maximal energy, which can be reached through 
astrophysical accelerating processes in a bottom-up scenario, and ^mm/max take into account that 
nearby /very early there are no astrophysical sources. As we will see in our comparison with 
UHECR data in Sect. EJ the overall normalization jo is fixed by the observed flux, and our 
predictions are quite insensitive to the specific choice for -E max , ^min, and z max , within their 
anticipated values. The main sensitivity arises from the spectral parameters a and m, for which 
we determine the 1- and 2-sigma confidence regions in Sect. EJ 

The injected protons propagate through the CMB. This propagation can be described [27] by 
Pb\ a ( r , E) functions, which give the expected number of particles of type b above the threshold 

2 This feature is shared with all strongly interacting neutrino scenarios. In contrast to these other scenarios, in 
our case, however, the threshold energy is automatically fixed by Standard Model parameters (Mw and aw)- 
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energy E if one particle of type a started at a distance r with energy With the help of 
these propagation functions, the differential flux of protons (b = p) and cosmogenic neutrinos 
{b = Vi,Vi) at earth, i.e. their number Nf, arriving at earth with energy E per units of energy, 
area (A), time (t) and solid angle (fl), can be expressed as 

In our analysis we go, according to dz = — (1 + z) H(z) dr/c, out to distances -R max corresponding 
to z max = 2 (cf. Ref. [2B])> while we choose z min = 0.012 in order to take into account the fact 
that within 50 Mpc there are no astrophysical sources of UHECRs. We use the expression 

n M (i + zf + n A ] (3) 



H\z) = Hi 



for the relation of the Hubble expansion rate at redshift z to the present one. Uncertainties of 
the latter, H Q = h 100 km/s/Mpc, with h = (0.71 ± 0.07) x^ij |2S|, are included. In Eq. ©, 
Qm and with fl M + Q A = 1, are the present matter and vacuum energy densities in terms of 
the critical density. As default values we choose Om = 0.3 and fl\ = 0.7, as favored today. Our 
results turn out to be pretty insensitive to the precise values of the cosmological parameters. 

We calculated P&| (r, E) in two steps, i) First, the SOPHIA Monte-Carlo program |3U] was 
used for the simulation of photohadronic processes of protons with the CMB photons. For e + e~ 
pair production we used the continuous energy loss approximation, since the inelasticity is very 
small (~ 10~ 3 ). We calculated the Pb\ a functions for "infinitesimal" steps (1 -j- 10 kpc) as a 
function of the redshift z. ii) We multiplied the corresponding infinitesimal probabilities starting 
at a distance r(z) down to earth with z — 0. The details of the calculation of the Pb\ a (r, Ei] E) 
functions for protons, neutrinos, charged leptons, and photons will be published elsewhere pTTj . 

The determination of the propagation functions took approximately one day on an average per- 
sonal computer. In this connection, the advantage of the formulation of the spectra (0) in terms 
of the propagation functions becomes evident. The latter have to be determined only once and for 
all. Without the use of the propagation functions, one would have to perform a simulation for any 
variation of the input spectrum (at, m, . . .), which requires excessive computer power. Since the 
propagation functions are of universal usage, we decided to make the latest versions of —dPb\ a /dE 
available for the public via the World- Wide- Web URL www.desy.de/~uhecr . 

As a check on our propagation functions, we have compared our predictions for the spectra (J2J) 
with the ones presented in Ref. for some specific values of the spectral parameters (a, m, 
E max , z min , z max , jo) and found quite good agreement. 



3 Spectrum of instanton-induced air showers 

In this section, we exploit a recent prediction of the electroweak instanton-induced parton-parton 
cross-section and determine the spectrum of instanton-induced air showers, which are initiated 
by the cosmogenic neutrino flux (J2J) impinging on the earth's atmosphere. 

Let us start with a review of the current knowledge about electroweak instantons. In the Standard 
Model of electroweak interactions (Quantum Flavor Dynamics (QFD)) there are certain processes 
which fundamentally can not be described by ordinary perturbation theory. These processes are 
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Figure 1: The electroweak instanton-induced neutrino-quark cross-section as predicted in 
Ref. j2E] for the whole range of the CM energy vl (left) and near the maximum (right). 

associated with axial anomalies and manifest themselves as anomalous violation of baryon plus 
lepton number (B + L) |32]- They are induced by topological fluctuations of the non-Abelian 
gauge fields, notably by instantons |33]. In Minkowski space-time, instantons describe tunneling 
transitions between degenerate, topologically inequivalent vacua. The corresponding tunneling 
barrier is given by the energy of the sphaleron [31] , an unstable static solution of the Yang- Mills 
equations, and of order M sp « nM w / aw ~ 10 TeV. The corresponding processes violate B + L 
according to the selection rule AB = AL = —3. 

It is generally accepted that such topological fluctuations and the associated B + L violating 
processes are very important at high temperatures [33] and have therefore a crucial impact on 
the evolution of the baryon and lepton asymmetries of the universe 3 . It is, however, still debated 
whether manifestations of such fluctuations - involving notably the associated production of 
0(1/ aw) ~ 30 W/Z-bosons in addition to the anomalously produced quarks and leptons - may 
be directly observed in high-energy scattering processes [21] . Despite considerable theoretical |H] 
and phenomenological (221112] efforts, the actual size of the cross-sections in the relevant, tens of 
TeV energy regime was never unanimously established (for reviews, see Refs. |37| I43|). 

There is a close analogy [Hj between QFD and hard QCD instanton-induced processes in deep- 
inelastic scattering [13] . Recent information about the latter - both from lattice simulations jlH] 
and from the HI experiment at HERA [17] - has been used by one of the authors to learn about 
the fate of electroweak B + L violation and associated multi-W/Z production at high energies [2E] 
(for a review, see Ref. [IE])- The prediction for the electroweak instanton-induced neutrino-quark 
cross-section is displayed in Fig. Q]as a function of the neutrino-quark center-of-mass (CM) 

energy y/S. At small CM energies, the cross-section is really tiny, e.g. a^/J ~ 10~ 141 pb at 
\f§ ~ 3 TeV, but steeply growing. Nevertheless, it stays unobservably small, a^) < io~ 26 pb for 
\/§ < 22.5 TeV, in the quite conservative fiducial kinematical region inferred via the QFD-QCD 
analogy from lattice data and HERA. It was noted that a slight extrapolation towards larger 
energies - still compatible with lattice results and HERA - points to a cross-section 1CT 6 pb 
at a CM energy of about 30 TeV, which is within the reach of the Very Large Hadron Collider. 

3 Standard Model electroweak baryogenesis seems excluded, however, due to the weakness of the electroweak 
phase transition (for reviews, see Refs. [371|3H])j while thermal leptogenesis 13^100] is quite successful. 
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Figure 2: Left: Prediction of the electroweak instanton- induced neutrino- nucleon cross-section 
crf^v (solid) in comparison with the charged current cross-section a u c N (dotted) from Ref. [SU], as 
a function of the neutrino energy E v in the nucleon's rest frame. Right: Neutrino interaction 
length due to combined effects of charged current interactions and instanton-induced processes. 

In this Letter, we will use the prediction from Ref. even at higher energies 4 , up to and above 
\/| m 40 TeV, where the cross-section reaches its maximum of order a few millibarn (cf. Fig. ^ 
(right)). The corresponding neutrino-nucleon cross-section is obtained after folding the parton 
cross-section a^y with the quark density functions f q , 



where s denotes the neutrino-nucleon CM energy squared and \i the factorization scale. For our 
numerical integration we have used various sets of parton distributions as they are implemented 
in the parton distribution library PDFLIB jl^j. Uncertainties associated with different parton 
distribution sets are in the O(20) % range and are not explicitely shown in the following. Figure 121 
(left) displays the prediction of the electroweak instanton-induced neutrino-nucleon cross-section 
as a function of the neutrino energy E v in the nucleon's rest frame for a choice \x = M\y of the 
factorization scale. Above a threshold at about E v ~ 10 18 eV, it quickly reaches one millibarn at 
about 10 19 eV, and tends to grow power-like, due to the growth of the sea quark distributions in the 
nucleon at small x, quite analogous to the standard charged current cross-section a™ N (cf. Fig. |2] 
(left)). It should be noted that such a cross-section will lead, via dispersion relations, to lower 
energy deviations of Standard Model predictions for elastic scattering from their perturbative 
values j2Zj- However, it is easily checked that, for the one shown in Fig. El (left) , these corrections 
will be unobservably small in the energy regime available at present accelerators 

The corresponding neutrino interaction length \ v = rrip/al ^, with = af N + crfy, is shown in 
Fig. El (right). It falls below Xo = 1031 g/cm 2 - the vertical depth of the atmosphere at sea level 5 - 
for E v > 3 - 10 19 eV. The apparent success of our scenario is based on the unexpected coincidence of 

4 It should be kept in mind, however, that, at the energies of interest here, the prediction in Fig. ^ is rather an 
educated extrapolation or guess (cf. Ref. |2S])- 

5 For our numerical calculations involving the atmospheric depth X(9) we have used a parametrization of the 
US Standard Atmosphere (1976) from Ref. [S3 - 
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this scale and E GZK . Above this energy, the atmosphere becomes opaque to cosmogenic neutrinos 
and all of them will end up as air showers. Quantitatively, this fact can be described by 



I — e 



(5) 



which gives the spectrum of neutrino-initiated instanton-induced air showers, for an incident 
cosmogenic neutrino flux F u = Y^i[F Ui + F^] from Eq. (J2J), in terms of the atmospheric deptff^ 
X(8), with 9 being the zenith angle. For E u >3 ■ 10 19 eV, one has \ V (E V ) < Xq, and the 
spectrum © quickly equals the incident cosmogenic neutrino flux, F V (E). For E u <4 ■ 10 18 eV, 
on the other hand, the cross-section cr**^ (.£/„) ^ 0.56 mb corresponds to a neutrino interaction 
length \ v {Ep) > 3000 g/ cm 2 , which is comparable to the atmospheric depth at larger zenith angles, 
8>70°. Therefore, for these energies, neutrino- initiated electroweak instanton-induced showers 
can be searched for at cosmic ray facilities by looking for quasi-horizontal air showers, 9 > 70° |25j . 
At the end of Sect. EJ we will show that the rate from our prediction (0) is consistent with 
observational constraints found by the Fly's Eye jHSl and AGASA (33] collaborations. 

Equation (j3J) does not account for the efficiency of an air shower array to trigger on low altitude 
air showers. Below 10 19 eV, neutrino-induced showers may be initiated so close to the array that 
the showers do not spread out sufficiently to trigger the array. As discussed in Ref. [25] , one 
may suppose that an array does not trigger on showers initiated within X tT = 500 g/cm 2 of the 
detection level. This can be implemented in Eq. (jSJ) by replacing X(6) with (X(6) — X tr ). Such 
an assumption seems reasonable for vertical showers seen by a ground array (AGASA), but is 
somewhat pessimistic for showers at larger zenith angles or for fluorescence detectors (HiRes). 
We have performed our fit in Sect. |U with/without such a "trigger" cut for AGASA/HiRes data. 
Its effect, however, turned out to be negligible. 

Proton-initiated electroweak instanton-induced air showers have been quite intensively studied 
in Ref. [23] and compared to generic proton- or iron- initiated air showers 6 . While identifiable 
systematic differences between average showers of different origin could be found, the differences 
did not appear to be sufficient to discriminate between proton-initiated instanton-induced showers 
and fluctuations in generic showers. The same is expected for neutrino-initiated instanton-induced 
air showers, as long as the first interaction occurs sufficiently high in the atmosphere, at a depth 
<500 g/cm 2 , which happens in our case for E^^IO 20 eV. We will find from our fits in the 
next section that the contribution of cosmogenic neutrino-initiated air showers to the UHECR 
spectrum starts to dominate at around this energy over the proton-initiated generic component. 



4 Comparison with UHECR data 

In this section, we compare the predicted air shower spectrum from Eqs. (J2J) and (0), the latter 
averaged over the appropriate range of zenith angles 9, 

Fpred{E] Of, 771, E TaaX) Z m i n , £ m ax; jo 

)=F p (E;...) + F^(E;...) } (6) 

6 On account of <§; cr^v w 100 mb, where <7p™ is the cross-section for generic proton-nucleon processes, the 
contribution of proton-initiated instanton-induced air showers to the UHECR spectrum can be safely ignored. 
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Figure 3: Ultrahigh energy cosmic ray data (points with statistical error bars; left: combination 
of Fly's Eye and HiRes data; right: combination of Akeno and AGASA data) and their best fits 
within the electroweak instanton scenario (solid) for E max = 3 ■ 10 22 eV, z min = 0.012, z max = 
2, consisting of a proton component (dotted) plus a cosmogenic neutrino-initiated electroweak 
instanton-induced component (dashed). 



with the observations. We perform a detailed statistical analysis and present a measure for the 
goodness of the instanton-induced scenario. 

The analysis consists of two parts, i) The UHECR collaborations give their results for the 
incoming flux in a binned form. Note, however, that the number of events in a given bin is integer 
and follows the Poisson distribution. In order to be able to give the goodness of the instanton- 
induced scenario by statistical methods, we determine the number of experimentally observed 
events in a given energy bin by converting the published values of the cosmic ray flux. We analyse 
the results from different experimental settings separately and perform the UHECR analysis for 
the two most recent results from the HiRes and AGASA collaborations, respectively, ii) We 
determine the 1-sigma and 2-sigma confidence regions for the parameters (a,m) characterizing 
the proton injection spectrum (cf. Sect.|2J). The method is similar to the frequentist's analysis [2H] 
and uses a Monte-Carlo integration in the mult i- dimensional space of bins. 

ad i) In our comparison, we use the observed data from \og(E/eV) = 17.2 to \og(E/eV) = 21. We 
have altogether 38 bins. The bins with the largest energies are empty. This non-trivial information 
is incorporated into the analysis, too. In the low energy region, there are no published results 
available from AGASA and only low statistics results from HiRes-2. Therefore, we included 
the results of the predecessor collaborations - Akeno [53] and Fly's Eye, respectively - into 
the analysis. With a small normalization correction, it was possible to continuously connect the 
AGASA data pQ with the Akeno ones and the HiRes- 1 monocular data f2| with the Fly's Eye stereo 
ones |2j, respectively (cf. Fig. EJ). Usually, it is advisable to avoid the combination of different 
experimental data. Since in the present case it is interesting to see how well our scenario works 
for energies below and above the GZK cutoff, we used the less problematic solution and combined 
results from experiments with the same techniques and with largely overlapping experimental 
groups. The normalization was matched at log(E/eV) = 18.5 for both cases. 

Note that the highest energy event of HiRes was published using a five times bigger bin size than 
for other energies 0. In order to preserve information, we prefer to keep the binning, give the 
particular bin with one event, and present upper bounds for the bins with zero event. From the 
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published data of HiRes, we determined the approximate energy and the corresponding bin of 
the highest energy event. 

ad ii) The goodness of the scenario is determined by a statistical analysis. In order to give the 
confidence region in the a-m plane, we determined the compatibility of different (a, m) pairs with 
the experimental data. For some fixed (a,m) pair, one can determine the expected number of event 
in individual bins (A = {Ai, A r }, where the Aj-s are non-negative, usually non-integer numbers 
and in our case the number of variables corresponds to the number of bins, thus r=38). For this 
specific (a, m), the probability distribution in the i-th bin is given by the Poisson distribution with 
mean Aj. The r dimensional probability distribution -P(k) is just the product of the individual 
Poisson distributions (here k = {k\, ...k r } is a set of non-negative integer numbers). It is easy 
to include also the ~ 30% overall uncertainty in the energy measurement of the experiments 
into the P(k) probability. We denote the experimental result by n = {ni, ...n r }, where the n^-s 
are non-negative, integer numbers. According to the r dimensional probability distribution, the 
experimentally observed event set n has a definite, though usually very small probability -P(n). 
The (a, m) pair is compatible with the experimental results if 

E m <s. (7) 

k|P(k)>P(n) 

For a l-(or 2-)sigma compatibility one takes s=0.68 (or s=0.95), respectively. The best fit is found 
by minimizing the sum on the left hand side. This technique is equivalent with the x 2 technique 
for a large class of problems 7 . Note, however, that the y 2 technique always gives a confidence 
region and the xVd.o.f is used as an estimate for the goodness of the scenario. Since xVd.o.f 
can be directly interpreted for Gaussian problems only, our goodness of the scenario technique is 
more general. 

Since we have 38 variables, it is practically impossible to calculate the sum in equation (J7J) exactly. 
Fortunately, there is no need for the exact calculation, the sum can be determined with arbitrary 
precision by using an importance sampling based Monte-Carlo summation. Since the sum of the 
individual probabilities is one, the left-hand-side of equation ((7J) can be rewritten as 

p(k) = ^- p(q) 9[p{k) ~ p(p)] 



k|P(k)>P(n) 



Ek^(k) 



Equation (jHJ) defines the Monte-Carlo summation straightforwardly. When calculating the sum, 
numbers with Poisson distribution are generated for k and only those are taken in the sum, for 
which P(k) > P(p). 

Figure H3 shows our best fits for the HiRes and for the AGASA UHECR data (the lower energy 
data is also included, as we explained before). The best fit values are a = 2.68, m = —0.1, for 
HiRes, and a = 2.68, m = —0.35, for AGASA. We can see very nice agreement with the data 
within an energy range of nearly four orders of magnitude. The fits are insensitive to the value 
of E max as far as we choose a value above ~ 3 • 10 21 eV. The shape of the curve between 10 17 eV 
and 10 19 eV is mainly determined by the redshift evolution index m. At z — 0, below 10 18 eV the 
attenuation length of protons is already around the size of the universe. Therefore, one would 
expect no distortions of the injected spectrum below this energy and an accumulation of particles 



7 For the application of the x 2 technique with UHECR data, see Refs. [Ill . 
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Figure 4: Confidence regions in the power-law index a and the redshift evolution index m of 
the primary proton injection spectrum, for fits to the Fly's Eye + HiRes data (1-sigma (solid); 
2-sigma (short-dashed)) and to Akeno + AGASA data (2-sigma (long dashed)), respectively, for 

-^max^3 ■ 10 eV, Z m [ n > 0, Z ma x = 2. 

just above it. However, at larger redshifts, the interaction lengths are smaller and the spectrum 
of particles created at cosmological distances has an accumulation peak at lower energies. The 
more particles are created at large distances (i.e. the larger m is), the stronger this effect will 
be 8 . The peak around 4 • 10 19 eV shows the accumulation of particles due to the GZK effect. 
Neutrinos start to dominate over protons at around 10 20 eV. 

It is important to note that, if we omit the neutrino component, then the model is ruled out on 
the 3-sigma level for both experiments. This is due to the fact that there are no nearby sources 
(-2mm 7^ 0) and all the events above 10 20 eV are highly inconsistent with the predictions. A choice 
of Zmin = makes the HiRes data compatible with a proton-only scenario on the 2-sigma level 
(see also Refs. jH |H])- If neutrinos are included, then - as they dominate over protons above 
10 20 eV - the fit results cease to be sensitive to the value of 2 mm . 

Figure |U displays the confidence regions in the a — m plane for HiRes and AGASA. The scenario 
is consistent on the 2-sigma level with both experiments. For HiRes, the compatibility is even 
true on the 1-sigma level. It is important to note that both experiments favor the same values for 
a and m, demonstrating their mutual compatibility on the 2-sigma level. If we ignore the energy 
uncertainty in the determination of the goodness of the fit, they turn out to be inconsistent. 

Finally, let us discuss the consistency of our scenario with the currently available limits on deeply 
penetrating showers from Fly's Eye and AGASA [S3]. Taking into account - in distinction 
to Ref. [HZ] - the atmospheric attenuation of the cosmogenic neutrino flux predicted in our 
scenario and the uncertainties in the estimate of the range of depth within which the shower 
must originate to trigger the array, we find that AGASA should have seen 1 10 quasi- horizontal 
air showers (#>60°) from the electroweak instanton-induced processes during a running time of 
1710.5 days. This is consistent with AGASA's present analysis of their respective data |53| . The 
Fly's Eye upper limit on the product of the total neutrino flux times neutrino-nucleon cross- 
section, {F v cr^jy) Fiy's Eye EH] j m the energy range 10 17 ^ 20 eV, can be translated, for a given 

8 Our finding suggests that the extragalactic UHECR component begins to dominate over the galactic one 
already at » 10 17 eV. If we start our fit at 10 185 eV - assuming that the galactic component dominates up to this 
energy - we find a very mild dependence on m and the same best fit values for a, with a bit larger uncertainties. 
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predicted neutrino flux FP red , into an upper limit on < {F v cr*^) Fly / s Eye /FP red , as long as it is 
smaller than 10 fib l2*ol I58j. We find that, for our predicted cosmogenic neutrino flux, the right- 
hand-side of this inequality is larger than 10 /u,b in the whole energy range, such that the Fly's 
Eye non-observation of quasi- horizontal air showers does not give any constraint. We therefore 
conclude that our prediction of the neutrino-nucleon cross-section, as shown in Fig. El (left), does 
not contradict any constraints from cosmic ray experiments so far, as long as the ultrahigh energy 
cosmic neutrino flux is at the cosmogenic level we have predicted. 

5 Summary and conclusions 

We have shown that a simple scenario with a single power-like injection spectrum of protons can 
describe all the features of the UHECR spectrum in the energy range 10 17 ^ 21 eV. In this scenario, 
the injected protons produce neutrinos during their propagation by interacting with the CMB. 
Through Standard Model electroweak instanton-induced processes, these neutrinos may interact 
with the atmosphere and give rise to a non-negligible contribution to the detected air showers at 
the highest energies. The model has few parameters from which only two - the power index a and 
the redshift evolution index m - has a strong effect on the final shape of the spectrum. We found 
that for certain values of a and m this scenario is compatible with the available observational data 
from the HiRes and AGASA experiments (combined with their predecessor experiments, Fly's 
Eye and Akeno, respectively) on the 2-sigma level (also 1-sigma for HiRes). The ultrahigh energy 
neutrino component can be experimentally tested by studying the zenith angle dependence of the 
events in the range \0 l8 + 20 eV at cosmic ray facilities such as the Pierre Auger Observatory and by 
looking for spatially compact energetic bundles at neutrino telescopes such as AMANDA [53] . 

Finally, let us emphasize that the same fit results are valid for all strongly interacting neutrino 
models if the neutrino-nucleon cross-section has a similar threshold-like behaviour as in Fig. |2l 
The instanton scenario, however, has the advantage that it is based solely on the Standard Model 
and can be falsified in the near future by a negative outcome of QCD instanton searches at HERA. 
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